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We present the optimal phase estimation for qubits in mixed states, for an arbitrary number of 
qubits prepared in the same state. 
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I. INTRODUCTION 

The problem of measuring the quantum phase has been 
a very long standing one in Quantum Mechanics, since 
London's[lj and Dirac's first attempts [2j in the late twen- 
ties. One of the main motivations is that the estimation 
of the phase shift experienced by a quantum system is 
the only way of measuring time with high precision in 
Quantum Mechanics, since we lack a time observable. 
This posed the problem of quantum phase estimation 
naturally within the framework of frequency standards 
based on atomic clocks ; 3|, and more generally, in high 
precision measurements and interferometry, the typical 
scenario in which the sensitivity of phase estimation is 
profitably used. 

More recently, the encoding of information into the 
relative phase of quantum systems is exploited in quan- 
tum computation and communication. In fact, in quan- 
tum computing most of the existing quantum algorithms 
can be regarded as multiparticle interferometers, with the 
output of the computation encoded in the relative phase 
between different paths [4|]. On the other hand, in some 
cryptographic communication protocols (e. g. BB840) 
information is encoded into phase properties. 

The above numerous applications had focused a great 
deal of interest on the problem of optimal phase esti- 
mation, which has been widelystudied in a thousand of 
papers (see for example Ref. || ) since the beginning of 
quantum theory [ij, |2( • The first satisfactory partial so- 
lution of the problem appeared in the late 70's (see Refs. 
and @ for reviews), and these works are generally 
regarded as one of the major successes of quantum es- 
timation theory and covariant measurements, allowing a 
first consistent definition of phase, without the problems 
suffered by the original definition proposed by Dirac [2j in 
terms of an alleged observable conjugated to the number 
operator of the harmonic oscillator. 

In the covariant treatment of Ref. Q the estimated 
parameter is a phase shift resulting from the action of a 
circle group of unitary transformations, with generator 



a selfadjoint operator with purely integer spectrum. A 
generalization of this method to degenerate phase-shift 
generator has been presented in Ref. 0- Such general 
approach can be applied to any input pure state, along 
with a restricted class of mixed states, the so called phase- 
pure states P. Hoj|. 

The possibility of efficiently estimating the phase for 
mixed states is of fundamental interest for practical im- 
plementations, in the presence of unavoidable noise which 
generally turns pure states into mixed, and for estima- 
tion of local phase-shift on entangled states. As a matter 
of fact, the freedom in the choice of the optimal mea- 
surement which results from degenerate shift operators 
opens the problem of the stability of the quality of the es- 
timation with increasing mixing of the shifted state. The 
problem of optimal phase estimation on mixed states is 
also very relevant conceptually, the phase being one of the 
most elusive quantum concepts. The main reason why 
the problem of optimal phase estimation on mixed states 
has never been addressed systematically so far is due 
to the intrinsic technical difficulties faced in any quan- 
tum estimation problem with mixed states. In this paper 
we derive the optimal measurement for phase estimation 
on qubits in mixed states, for an arbitrary number N 
of qubits prepared in the same state, using either the 
Uhlman fidelity or the periodicized variance as a figure 
of merit. 



II. THEORETICAL DERIVATION 

Let us consider a system of N identical qubits prepared 
in the same mixed state pfi = ^(1 + n ■ a), where \n\ = 
r < 1 and oi are the three Pauli matrices. The total 
state of the N qubits is described by the density matrix 
Rft = P® N ■ The phase transformation U$ is generated 
by the z component of the total angular momentum J z — 
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ai^ , namely 
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The problem is now to estimate the unknown phase-shift 
</> on the known state Rf t . We consider a covariant mea- 
surement, namely we require that the efficiency of the 
measurement procedure does not depend on the value of 
the phase to be estimated. In this case, without loss of 
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generality we can assume that the initial state pa has 
no component along a y , corresponding to real matrix p^ 
in the a z representation. The phase estimation problem 
then resorts to find the best POVM P(d </>) for deter- 
mining the unknown parameter (j) in Eq. Q). The fact 
that P(dcj)) is a POVM corresponds to the constraints 
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In the quantum estimation approach the optimality is 
defined by maximizing the average of a given figure of 
merit C (</>,</>'), assuming a uniform prior distribution of 
the parameter <fi 



(C) = 
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C&MTrpiRaUlPid </>')], (3) 



where C(4>, <f>') = C(<f>— (/>'). In Ref. y| it was proved that 
the solution for an estimation problem covariant under a 
unitary group representation can be written as the group 
orbit under the same representation of a fixed positive 
operator £ (called seed of the POVM), and for the present 
case one has 



P(d< 



(4) 



In the following we will denote by \m, a) an orthonormal 
basis, with to denoting the eigenvalues of ^ J z , which la- 
bel the equivalence classes of irreducible representations 
of the group {U^}, while a is a degeneration index, cor- 
responding to the multiplicity space of the representa- 
tion to. The normalization condition (J2J for the POVM 
P{d<p) implies that (m,a\£\m,b) — for a ^ b, whereas 
(to, o|£|to, a) = 1 for all a. 

In quantum estimation theory the quantity to be min- 
imized can be always written in the form of the expecta- 
tion of a cost operator as follows 



(5) 



The choice of the cost function C(<p) depends on the es- 
timation criterion. The most commonly adopted criteria 
are the periodicized variance 



C(4>) = v(</>) = 4 sin 2 | = 2(1 - cos ( 



(6) 



or the (opposite of the) fidelity between the true and 
the estimated states, which for mixed states has the well 
known Uhlman's form 111] 



C{<t>) = l-F{4>) = l- 



Tr 



which for qubits simplifies as follows |l^ 



l-P(^) = y(l-cos< 
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Both cost functions depend on (f> only through its cosine, 
whence we need to maximize the averaged cos <j), namely 
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(j*' +e- i *)Tr[R n U^Ul,]. (9) 



The evaluation of the integral in Eq. © leads to the 
following expression 

(c) = Re ^ (m,a\£\m+ 1, b)(m + 1, b\R H \m, a) . (10) 



We now decompose p® N into irreducible representations 
of §U(2), as shown in Ref. 0] recasting Rft into invariant 
block-diagonal form on the orthonormal basis \ j, m,a)r — 
Uj i0l \j,m, l)g for the minimal invariant subspaces of the 

§U(2) representations, with b = — , and Uj, a denoting a 
suitable set of unitary operators. 

Rh=pT= E M'E^ii- (H) 

j=((iV/2» a=l 



m——j 

\j,m,l) t = |j,m) s ® \*-)® J -\ 
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where r± = j(l ± r), ((x)) denotes the fractional part of 
x (i. e. ({N/2)) = for N even and {{N/2)) = 1/2 for 
N odd), J = N/2, and dj is the multiplicity of the j-th 
irreducible representation of SU(2) 
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whereas I*-) denotes the singlet state. This decomposi- 
tion is useful since to, j, a label also the irreducible repre- 
sentations of {Uff,}, to being the eigenvalue of J z , and j, a 
becoming both degeneration indices. The block diagonal 
form of Rft shows that the only coupling produced by 
the phase shift between irreducible components with to 
and m + 1 can occur only between vectors in the same 
invariant subspace j,a of SU(2). Upon recasting Ra in 
the form of Eq. pi|) , the value of (c) in Eq. i|lU|l involves 
only the following terms 

(c) = Re ^2 (m,ja\£,\m + l,ja)(m + 1, ja\R n \m, ja) , 



(15) 

where we used the short notation \m,ja) = \j,m,a) z , 
since the subspaces j, a are invariant under any unitary 
in SU(2), and \j,m, a)g = T^\g)\j,m,a) z for some g € 



§U(2). Now, the following bounding hold 



(c) < }^ (m,j, a\(\m + l,ja){m + 1, ja\Rfj\m, ja) 

7i, j a 

< X] K m ' j' a l£l m + l,ja\Rn\m,ja)\ 

m,ja. 

< ^2 \(m+l,ja\Rfi\m,ja}\ , (16) 



where the last bound follows from positivity of £. We 
show now that all bounds can be achieved by a suitable 
choice of the operator £ compatible with constraints J2J- 
The first two bounds can indeed be achieved by choosing 
the phases of the matrix elements (to, j, a\£\m + l,ja) 
in such a way that they compensate the correspond- 
ing phases of (to + l,ja\Rft\m,ja). The last bound is 
achieved by just taking the moduli of the matrix elements 
(m, ja\£\m+l, ja) to be 1. It remains to prove that these 
choices are compatible with positivity. In order to show 
this, let us write 



(m+l,ja\Rfi\m,ja) = \(m+l, ja\Rft\m, j a >|e i *< m+1 ' mj ' Q > 

Since only the elements on the first over-diagonal and 
under-diagonal are involved, one can write the phases 
x(m + 1, m,ja) as the difference of two functions as fol- 
lows 



X(m + l,m,ja) = 7(771, ja) - 7(771 + 1, ja) , (18) 



as the number of independent linear equations in Eq. 
(|18|l is 2 N — 1 while the unknown phases are 2 N . Then 
one can take 
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mixed initial states, we have 
(c) = X \(m + l,ja\Rft\m,ja) 

= X ( r + r -) J x 



£ E ^+r_) J x 
j={{N/2)) m=-j 
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Notice that since we assumed that ri has no component 
along the direction y, then is just the rotation around 

the axis y connecting the oriented z axis with b, namely 

T(g^) = e i9J y for some 6, with J y = \ Y2 =1 <4*° • 



III. NUMERICAL RESULTS 

The expression for the Wigner matrix elements 
T if) (gz)ik is given by 



(j + Z — -fc-i)!(t-Z + fc)!t! 

^sin 2t -'+^. 
2 2 

(22) 



cos 2 ^+ i - fe - sin 2t -'+ fe 



The explicit expression of Eq. (|21() is very lengthy, and 





where |e(j, a)) is the generalized Susskind-Glogower vec- 
tor 

\e(j,a))= J2 e^ m ^\m,j,a). (20) 

m=-j 

It is immediate that Eq. (|19|) represents a positive opera- 
tor and by construction £ produces a normalized POVM, 
while achieving the bounding in Eq. (|16[) 

Specifically, for a collection of identically prepared 



FIG. 1: The plot of (c) as a function of 8 and r. The plots 
correspond to systems of 10 qubits and 20 qubits respectively. 

has been evaluated using symbolic calculus for J up to 
21/2, namely for a total number of spins equal to 21. 
The plot of the averaged cosine (c) as a function of 9 and 
r is represented in Fig. ^ an d exhibits two interesting 
intuitive features. The first is that the maximum versus 
9 occurs for 9 = |, namely for qubits lying in the equa- 
torial plane. The second is the improving figure of merit 
versus the purity r. Equatorial pure qubits are optimal 
for phase detection, however, the figure of merit is quite 
stable around its maxima, still with N — 10 copies. 
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FIG. 2: The logarithmic plots represent 2(1 — (c)), where (c) 
is the averaged cosine, as a function of the number of spins 
N, for 0=1 and for the following values of r: ♦ r = .7, 
★ r = .8, ■ r = .9, A r = 1. 



[Ttij . there should exist a separable strategy (such as an 
adaptive scheme) which is not necessarily covariant, nev- 
ertheless it would be able to achieve the same Cramer- 
Rao bound asymptotically: such non covariant schemes, 
e. g. homodyne-based estimation of the phase, will be 
the subject of further studies. 
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Fig. [2]shows the averaged cosine (c) versus the number 
of qubits N for equatorial states. Numerically, for N — > 
oo we find the asymptotic behavior 2(1 — (c)) oc TV -1 . 
More precisely, for the Uhlman fidelity F in Eq. J7J we 
find an asymptotic behavior saturating the Cramer- Rao 
lower bound [15j. This gives a strict lower bound for 
variance Ac/) 2 valid for any estimate. For independent 
copies, one has 



1 



Acf 2 > —^[{dp/d^C]- 1 



N 



(23) 



where for each (j) the operator C is defined by the identity 



dp/d^ = -{ P c + Cp). 



(24) 



Notice that the bound holds for any estimate, whence re- 
gardless the nature of the measurement (corresponding 
to either joint or separable POVM's). Since the estima- 
tion is covariant, we can just consider (j) — 0. A simple 



evaluation shows that C 
then given by ^ Tr[p £ 2 ]" 

A0 2 > 



r cos 9a. 



Nr' z cos- 



and the bound is 
-s, namely 
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N r 2 cos 2 



(25) 



For small A(f> 2 using the Taylor expansion of the co- 
sine one has A0 2 ~ 2(1 — (c)). In Fig. [3] we plot 
2(1 — (c))N of our optimal estimation for 9 = ver- 
sus r for TV = 16, 18, 20, against the Cramer-Rao bound 
p-. From the comparison we see that our estimation 
approaches the Cramer-Rao bound for large N. Notice 
that according to recent studies of theoretical statistics 



FIG. 3: The logarithmic plot of 2N(l - (c» vs r, for N = 
16, 18, 20 and 8 = 0. The line on the bottom represents the 
bound given by the Cramer- Rao inequality, namely 1/r 2 



IV. CONCLUSIONS 

In conclusion, we have presented the optimal measure- 
ment for phase estimation on TV qubits all prepared in 
the same arbitrary mixed state. The Uhlman fidelity 
saturates the Cramer-Rao bound for this problem, con- 
firming the optimality of the measurement. The optimal 
estimation is achieved for equatorial qubits and generally 
the fidelity is improving with purity. The specific form of 
the optimal POVM in terms of the generalized Susskind- 
Glogower vector in Eqs. (|2Tj|) and lfH5|) suggests possible 
physical implementations in terms of a generalized mul- 
tipartite Bell measurement. 
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